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ABSTRACT 

We present a new planetary nebula, first identified in images from the Australia Telescope 
Compact Array, although not recognized at that time. Recent observations with the Spitzer 
Space Telescope during the GLIMPSE Legacy program have rediscovered the object. The high- 
resolution radio and infrared images enable the identification of the central star or its wind, the 
recognition of the radio emission as thermal, and the probable presence of polycylic aromatic 
hydrocarbons in and around the source. These lead to the conclusion that G313.3-1-00.3 is a 
planetary nebula. This object is of particular interest because it was discovered solely through 
radio and mid-infrared imaging, without any optical (or near-infrared) confirmation, and acts 
as a proof of concept for the discovery of many more highly extinguished planetary nebulae. 
G313.3-1-00.3 is well-resolved by both the instruments with which it was identified, and suffers 
extreme reddening due to its location in the Scutum-Crux spiral arm. 

Subject headings: planetary nebulae — radio continuum: general — infrared: general — stars: winds 


1. Introduction 

The GLIMPSE (Galactic Legacy Infrared Mid- 
Plane Survey Extraordinaire) Legacy program 
(Churchwell et al. 2004) is surveying 220 deg^ 
of the plane in all four bands of the Infrared Ar¬ 
ray Camera (IRAC: Fazio et al. 2004) aboard 
the Spitzer Space Telescope (Werner et al. 2004). 
GLIMPSE is generating a wealth of detailed im¬ 
agery of the plane at mid-infrared (MIR) wave¬ 
lengths. The greatest scientific yield from any 
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such survey occurs when multiwavelength data 
are compared. Detailed panoramic comparisons 
of areas of the Galactic plane in the radio contin¬ 
uum and MIR were undertaken by Cohen & Green 
(2001), who examined 8 deg^ in the I = 312° re¬ 
gion using matched resolution imagery. These 
authors investigated the detailed 843-MHz inven¬ 
tory of this field made by Whiteoak et al. (1994) 
with the Molonglo Observatory Synthesis Tele¬ 
scope (MOST; Mills 1981; Robertson 1991). The 
Midcourse Space experiment (MSX: Price et al. 
2001) detected a substantial number of Hii re¬ 
gions, of diverse morphology, in its MIR survey of 
the Galactic Plane, which Cohen & Green (2001) 
used to calibrate the ratio of 8.3-p,m to 843-MHz 
flux densities, F8.3Aim/5'843M7rz, as a discriminant 
between thermal and nonthermal radio sources. 
We are engaged in a detailed comparison of the 
GLIMPSE images of the same I = 312° field pre¬ 
viously studied with MSX and the MOST. 

This paper presents high-resolution radio con- 
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tinuum and Spitzer images of G313.3+00.3, an 
object which appears as a nebulous ring with a 
central star. The radio images, together with a 
new 843-MHz image, have been used to estimate 
a spectral index for the source, which implies that 
the object is thermal (§2). Previous infrared (IR) 
and optical information available in the literature 
constitute §3. IR images of the nebula taken with 
the IRAC are discussed in §4, and alternative ex¬ 
planations for the nebula are considered in §5, af¬ 
ter which (§ 6 ) we conclude that the source is an 
extremely highly reddened planetary nebula (PN). 

2. Radio continuum observations 

The original impetus for carrying out high- 
resolution radio continuum images of this region 
was to search for posssible counterparts to an 
unidentified Galactic plane 7 -ray source. The ini¬ 
tial radio survey used to pinpoint the target region 
was the first epoch Molonglo Galactic Plane Sur¬ 
vey (MGPSl; Green et al. 1999) from the MOST. 
That 843-MHz continuum observation has a reso¬ 
lution of 43" X 43" cosecJ and a la sensitivity of 
1 — 2 mJy beam“^. The MOST image used in the 
present paper to calculate flux densities is a mosaic 
(J1430M60) of the second epoch survey MGPS2 
(Green 2002) which is available at the website 
http://www.astrop.physics.usyd.edu.au/SUMSS. 
The observational parameters are essentially un¬ 
changed from MGPSl. 

The targeted radio continuum images were part 
of the project to search for the radio counterpart 
of the EGRET source GeV J1417-6100 (Roberts 
et al. 1999). Observations were made at 13 and 
20 cm with the Australia Telescope Gompact Ar¬ 
ray (ATGA; Frater et al.l992) in three config¬ 
urations: 1.5A (1998 January 11), 0.75A (1998 
April 22), and 0.375 (1998 March 28). In ad¬ 
dition, high-resolution 20 -cm data were included 
from three observations taken in 1999 and 2001 in 
the 6 D, 0.375, and 0.75G configurations by Froney 
Grawford and Simon Johnston (these are detailed 
by Roberts et al. (2001)). The ATGA can ob¬ 
serve two simultaneous frequency bands and 1384 
and 2496 MHz were selected. For the 1.5A ar¬ 
ray, the lower frequency was slightly different, at 
1344 MHz. For all observations a bandwidth of 
128 MHz was used and full Stokes parameters were 
recorded. For most of the analysis in Roberts et 


al. (1999), images were used which did not include 
the baselines to the 6 -km antenna. 

High-resolution images including the 6 -km an¬ 
tenna of the ATGA were also made, using the 
MIRIAD (Sault et al. 1996), maximum entropy 
method to minimize artefacts. The ring source 
presented in this paper is shown in the high- 
resolution 13-cm image of Figure 3 of Roberts et 
al. (1999), but there is no discussion on its na¬ 
ture. This image has an angular resolution of 
2.8" X 2.5" and an rms noise of 0.1 mJy beam“^. 
The high-resolution 20-cm image presented in this 
paper has an angular resolution of 5.1" x 4.6" 
and an rms noise of 0.1 mJy beam“^. In both 
the high-resolution images a shell is clearly seen. 
In the low resolution images, the source is seen as 
slightly extended which allows an accurate mea¬ 
sure of the total flux density for spectral index cal¬ 
culations. These total flux densities at the three 
frequencies appear in Table 1. A formal regres¬ 
sion yields a slope of 0.000 ± 0 . 002 , and an essen¬ 
tially constant flux density of 25 ± 5 mJy, with 
= 4.3, consistent with a flat spectrum from 
thermal bremsstrahlung. Note that, at these fre¬ 
quencies, the Galaxy is optically thin so the is¬ 
sue of synchrotron self-absorption along the line- 
of-sight does not arise. 

Figure 1 presents the high-resolution 13-cm im¬ 
age of the nebula overlaid by contours of high- 
resolution 20-cm emission. A central point source 
is detected at both frequencies although its peak 
is better determined at 13 cm: 14'‘18”"27.489® ± 
0.007^ -60°47'10.49" ± 0.06" (J2000), corre¬ 
sponding to Z=313.3541°, &=-|-0.3125°. It lies at 
the geometric center of the nebula. We identify 
this as the location of the central star of the neb¬ 
ula. 

Figure 2 shows the spatial variation of spectral 
index between 13 and 20 cm across the nebula and 
at the stellar location. This map was made by 
matching the 13-cm image to the resolution and 
pixel grid of the 20-cm map. The ratio of these 
two images was scaled to represent the spectral 
index, a, where Sj^ oc At the star, a « 1.6, 
which is somewhat flatter than a Rayleigh-Jeans 
slope for emission from optically thick gas. At the 
location of the central star is a 13-cm point source 
with a flux density of 1 mJy. It is difficult to ex¬ 
tract a 20 -cm counterpart because locally there 
is extended emission of comparable brightness to 
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the stellar peak. The formal fit yields 2.6 mJy 
although this peak is ~1.5 beams away from the 
well-determined 13-cm peak. If this does represent 
the star at 20 cm then the spectral index is ~1.7, 
with probable uncertainty large enough to include 
an index of 2, i.e. emission in the Rayleigh-Jeans 
regime of a hot object. The region immediately 
surrounding the star is occupied by nebular gas 
with an average index of 0.75 ± 0.08. The onter 
annulus of the nebula has a = —0.09 ± 0.06, in¬ 
dicative of thermal (free-free) emission. 

3. Infrared and optical data from the lit¬ 
erature 

The first detection of a MIR source at this loca¬ 
tion was made by IRAS. The IRAS Point Source 
Catalog lists IRAS I4I47 — 6033, whose original 
(BI950) position corresponds to 14^18’"27.4^ ± 
2.6^ -60°47'08" ± 16" (J2000). It was detected 
only at 25 /rm with a flux density below 1 Jy. Con¬ 
sequently, the object has never been classified, be¬ 
cause neither color nor color-color criteria could 
be applied. 

During their MIR and radio study of the 
I = 312° region Cohen & Green (2001) used MSX 
imagery of the Galactic plane with 20" resolution 
from the “GB-02” data collection events. The 
nebula is faint but was detected in these images, 
although it was essentially unresolved. Additional, 
deeper MSX data are now available (the “CB-03” 
data) which are roughly 5 x deeper than the GB-02 
images. The nebula is not listed in the MSX Point 
Source Gatalog (Ver. 2.3: Egan et al. 2003), al¬ 
though it was covered in the GB-03 field centered 
near I = 315° (Price et al. 2001). In this deeper 
8.3-/im image it appears as a resolved object about 
30" X 30", i.e. with a deconvolved size of « 22". 
The same source is also detected in the 14.6-^m 
CB-03 image of this same field. These MIR flux 
densities are collected into Table 1. 

A search of the 2MASS archives and images 
shows no near-infrared (NIR) counterpart to the 
nebula. Likewise, no southern broadband opti¬ 
cal {Bj,R,I) image, from the UK Schmidt Tele- 
scope/SuperCosmos surveys shows a counterpart 
to the nebula. However, there may be an ex¬ 
tremely faint counterpart in the highly sensitive 
narrowband Southern Ha Survey of the Galactic 
plane (SHS: Parker & Phillipps 1998,2003), which 


can serve at least as an upper limit to Ha bright¬ 
ness (see §5). These attributes suggest a highly 
reddened and/or distant object. 

4. New Spitzer infrared images 

Cohen & Green (2001) established a distinc¬ 
tion between thermal and nonthermal radio emis¬ 
sion by comparing MIR and radio counterparts of 
many different types of Hii region in the I = 312° 
area. GLIMPSE offers a deeper exploration of 
this field and could potentially unravel the most 
complex areas in which MIR-emitting thermal fil¬ 
aments are juxtaposed with regions known to be 
nonthermal from their radio emission. When we 
examined the IRAC Band-4 (8.0-^m) mosaic im¬ 
ages in the vicinity of the EGRET source J1417 — 
6100 (Roberts et al. 1999) we immediately saw an 
IR ring nebula and a central star, where the lower- 
resolution MSX data showed just a faint blob at 
8.3 /rm. 

The individual IRAG frames were calibrated 
by the Spitzer Science Genter (SSG), and further 
processed by the GLIMPSE pipeline to remove 
instrumental artifacts such as column pulldown 
and banding (Hora et al. 2004a). Point sources 
were extracted from each frame using a modi¬ 
fied version of DAOPHOT (Stetson 1987), and 
cross-referenced using the SSC bandmerger. Mo¬ 
saics were created using the Montage package^. 
The final GLIMPSE image format has 0.6" pix¬ 
els and appears in Figure 3, which illustrates the 
four IRAC images of this nebula, which becomes 
progressively more conspicuous with increasing 
wavelength. The mosaic images conserve sur¬ 
face brightness. Spatially integrated, background- 
subtracted fluxes for the whole nebula were ex¬ 
tracted at each wavelength (see Table 1) after 
correcting for obvious point sources within the 
ring or lying near its rim. The flux calibration 
of GLIMPSE is cumulatively checked after com¬ 
pletion of each segment of the survey against a 
network of calibrators established in the survey 
area along the Galactic Plane (Gohen 2004). At 
present, the observed absolute flux densities of 
a sample of about 70 early A-dwarf and 60 K- 
giant calibrators show agreement with predictions 
to within ~3% in all IRAC bands (Cohen 2004). 


^available at http://montage.ipac.caltech.edu 
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At the center of the ring is a local peak in the 
IRAC Band 1 (3.6 /rm) and Band 2 (4.5 /im) 
images, which is swamped by bright nebulos¬ 
ity in Band 3 (5.8 /xm) and Band 4 (8.0 /xm). 
This peak pixel is at ^=313.3541°, &=-|-0.3125° 
and is clearly the IR counterpart of the central 
point source detected in the radio. Its desig¬ 
nation in the GLIMPSE archive is: SSTGLMA 
G3I3.3541-1-00.3124. Two circles with radii of 
12.6" and 7.6", centered on this NIR peak, pro¬ 
vide good fits to the outer and inner edges of the 
ring, validating the choice of the central source as 
the exciting star. PSF-fitting photometry of the 
central star was extracted from the 3.6 and 4.5- 
/xm images, above the locally determined nebular 
emission inside the ring. There is no noticeable 
central star in the 5.8 or 8.0-/xm images. 

Profiles in Galactic latitude through the cen¬ 
ter of the ring are presented in Figure 4. Each 
normalized profile was constructed from the av¬ 
erage of the three adjacent spatial slices centered 
on, and flanking, the position of the peak Bandl 
pixel of the central star. Negative distances from 
the center correspond to southward latitudes from 
the star. 

Figures 5 and 6 compare the high-resolution 
13-cm with the 4.5-/xm and 8.0-/xm images. In 
each IR band the lowest contour corresponds to 
4 ct above the local background level. Only the 4.5- 
/xm IRAC band should contain no significant con¬ 
tribution from PAH emission. Radiation in this 
band might represent line emission from ions (e.g. 
[Mgiv]), or molecular gas (e.g. ro-vibrational H 2 
lines, 0-0 S(8)-S(II)), thermal emission from free- 
free continuum, or from heated dust grains. 

The good match in nebular extent and mor¬ 
phology between the 4.5-yxm and radio maps would 
suggest that both these images trace ionized gas. 
By contrast, the 8.0-/xm emission clearly extends 
beyond the radio thermal emission, consistent 
with the presence of polycyclic aromatic hydro¬ 
carbon (PAH) emission in the 8.0-/xm bandpass. 
This would be expected if the PAHs resided in 
the photodissociation region (PDR) surrounding 
the ionized gas (e.g. Cohen & Green 2001). This 
8 .0-/xm extension beyond the ionized gas could 
also be consistent with the presence of molecular 
hydrogen emission lines in IRAC Band 4. 


5. The nature of the nebula 

We will calculate the mass and distance of the 
nebula, after determining the contribution of re¬ 
combination lines, and free-free and bound-free 
continua, to the IRAC fluxes, based on the ra¬ 
dio free-free fluxes, and estimating the interstellar 
extinction from an upper limit to the Ha flux, and 
from extinction maps based on DIRBE data. 

We have calculated the contributions suggested 
by the radio emission at 13 cm, assumed to be 
thermal, to the H-line emission within each of the 
IRAC filters. The strongest expected recombina¬ 
tion lines are Pfy and Pfd in IRACl; Bra and Pf/3 
in IRAC2; Huy in IRAC3; and Pfa and Hu/3 in 
IRAC4. All the 13-cm emission was attributed to 
free-free, using emissivities from Storey & Hum¬ 
mer (1995) for Te=10‘* K (the density is almost 
unimportant), approximating the radio free-free 
Gaunt factor by Osterbrock (1989: eqn. 4.30), 
and ignoring the effects of any dust mixed with 
the ionized gas. The free-free continuum within 
the IRAC bands was computed with IR Gaunt fac¬ 
tors from Beckert et al. (2000), and the bound-free 
component was estimated following Seaton (1960) 
for transitions from the continuum to levels with 
n = 6 to 9. The resulting intrinsic flux densi¬ 
ties expected from the entire nebula, combining 
these thermal processes, are 4, 8, 1, and 2 mJy in 
IRAC Bands 1,2,3,4, respectively. We now require 
an estimate of the intervening extinction to assess 
whether these processes make any observable MIR 
contributions to the nebular emission. 

The thermal radio emission provides an esti¬ 
mate of the intrinsic emission line flux in Ha (see 
Condon (1992), eqns. (3) and (4a)). Using this 
formulation yields 2.1x10“^^ erg cm“^ s“^ for the 
integrated intrinsic Ha flux, before the effects of 
extinction. To determine an upper limit for the 
spatially-integrated Ha flux of the nebula, we have 
applied the absolute calibration procedure to the 
SHS Ha image that Pierce (2004) describes. The 
SHS obtains short-red (SR) continuum exposures 
of each Ha field. Numerically scaled versions of 
these were intended to be used to remove any con¬ 
tinuum that is included in the narrowband Ha 
filter. After detailed cross-comparison with the 
absolutely calibrated, but spatially much lower- 
resolution, Southern Ha Sky Survey Atlas (Gaus- 
tad et al. 2001), Pierce provides both the SR scale 
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factors and the number of flat-fielded SHS Ha 
counts pixel” ^ that are equivalent to 1 Rayleigh 
(10”^ erg cm”^ s”^ sr”^ at Ha). The ring lies in 
the SHS field HAL0182, for which Pierce gives an 
SR scale factor of 1.12 and 7.8 counts pixel”^ 
although there is a very large amount of scatter 
about the latter value in this held (Pierce 2004, 
Table 2). Visual examination of this calibrated 
SHS Ha image, guided by overlaying the 13-cm 
outer radio contour of the nebula, suggests that a 
potential excess of Ha emission exists within the 
ring. However, there is no statistically significant 
excess compared with other regions of the image 
lacking point sources. Taking the measured Ha 
flux within the ring as a conservative upper limit 
to the Ha emission, the brightness integrated over 
the nebula is less than 9.0x10”^^ erg cm”^ s”^. 
Comparing this limit with the predicted Ha emis¬ 
sion implies that Ay exceeds 13™ (e.g. using the 
reddening law of Cardelli et al. (1989) in the op¬ 
tical, assuming Av/E{B — V) = 3.1 (appropriate 
for the diffuse interstellar medium)). Using the 
extinction maps of Schlegel et al. (1998), the to¬ 
tal Galactic extinction along the line-of-sight to 
the nebula is Ay~36™. Although this is large it 
is not atypical for the cores of spiral arms. 

The influence on IRAC flux densities of this 
substantial visual extinction cannot be neglected. 
Adopting an Ay of 36™ yields extinctions (±1 ct) 
in the IRAC bands of 2.28±0.25, 1.75±0.32, 
1.75±0.40, and 1.75±0.40™, respectively, from the 
empirical reddening law found by Indebetouw et 
al. (2005) for the IRAC bands along two lines-of- 
sight. A conventional reddening law (e.g. Cohen 
et al. 2003) gives extinctions of 1.85, 1.26, 0.94, 
and 1.59™, respectively. Both these sets of atten¬ 
uations would render the expected contributions 
of gaseous thermal processes to the observed flux 
densities negligible (1 mJy or less) in all IRAC 
bands. 

The radio flux density can be converted into 
the ionized mass of the ring (Mezger & Henderson 
1967, Appendix A, eqn. A.14), using Te = 10'^ K, 
assuming the simplest model of a uniformly filled 
sphere, and adopting a 5-GHz flux density of 
25 mJy, equal to the best fit to the three lower 
frequencies described above. The ionized mass is 
0.0022 Mq, where T)kpc is the distance to 

the ring in kpc. The mean ionized mass of plan¬ 
etary nebulae (PNe) is between 0.1 and 0.25 Mq 


(Boffi & Stanghellini 1994), suggesting that the 
ring is between 4.6 and 6.6 kpc from the Sun. In 
this direction the line-of-sight extinction based on 
reddened OBA stars amounts to only ~3™ out to 
a distance of 3 kpc (Neckel et al. 1980). The 
nebula is seen in the direction of the tangent to 
the Scutum-Crux spiral arm. The tangent point 
lies 6.6 kpc from the Sun (Taylor & Cordes 1993). 
We, therefore, favor a distance close to 6.6 kpc 
because it provides a natural explanation for the 
extremely high extinction that the PN suffers at a 
relatively modest distance from the Sun, namely 
in the dust lane associated with this arm. Identify¬ 
ing the reddening with an arm links the bulk of the 
DIRBE-derived extinction in this longitude to the 
spiral arm, rather than being distributed along an 
extended line-of-sight through the Galaxy as, for 
example, the 30™ of Ay associated with the direc¬ 
tion to the Galactic center. The implied nebular 
radius is 0.4 pc, close to the median (0.3 pc) for a 
set of PNe in the Magellanic Glouds (Stanghellini 
et al. 1999). 

6. Discussion 

The morphology of the MIR and radio contin¬ 
uum images suggests that ionized gas dominates 
the 4.5-/im image, while PAHs contribute dom¬ 
inantly at 8.0-/xm. However, as demonstrated 
above, the contributions to the 4.5-/im emis¬ 
sion from thermal processes that involve atomic 
hydrogen are very small. Therefore, any ionic 
lines within this band must come from heavy el¬ 
ements with rather strong gf values; for example, 
[Mgiv] at 4.485 /rm, and [Avi] at 4.530 /rm. Both 
these lines are seen in high-excitation PNe. In 
NGC 7027 and NGG 7662, which have significant 
MIR continua, the lines account for almost 10% 
of the flux in the 4.5-/im band. In a PN lack¬ 
ing such a continuum the lines’ contribution could 
be larger. Alternatively, a variety of H 2 emission 
lines could contribute in several IRAG bands, but 
chiefly at 8.0 /rm, as Hora et al. (2004b) found 
from an IRAC survey of nine PNe. 

Were the nebula to be a supernova remnant 
(SNR) we would expect a negative radio spec¬ 
tral index. Even if the SNR were young, such a 
shell source would have a negative index reflect¬ 
ing the synchrotron emission mechanism. Exam¬ 
ples of young SNRs with a ringlike structure in- 
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elude Tycho, Kepler, Cas A, and SN 1006, and 
all their radio indices are around —0.7. Accretion 
onto a compact object can produce a flat spec¬ 
trum, as could a plerion. The 13 — 20 cm spec¬ 
tral index (~1.7) for the central source is rather 
uncertain, but probably excludes a flat spectrum, 
making this interpretation unlikely. If the cen¬ 
tral source were a very young pulsar wind neb¬ 
ula, the spectrum could even be inverted as in 
SN 1986J, which recently showed a double power 
law with slopes —0.55 and -1-1.38 (Bietenholz et 
al. 2002). However, this time-dependent inversion 
appears only at high radio frequences (>10 GHz). 
These authors found synchrotron self-absorption 
within the source was an unlikely explanation for 
the steeply positive power law. The ring’s radio 
emission is too weak to measure polarization reli¬ 
ably with current instrumentation. 

The MIR morphology is also very different from 
that of a young SNR, e.g. Cas A. Figure 7 il¬ 
lustrates the MSX data for this young remnant, 
overlaying 21.3-/im contours on the 8.3-/im image. 
The SNR is highly clumpy, shows large appar¬ 
ent density contrast and ragged inner and outer 
edges. These characteristics are unlike those of 
G313.3-1-00.3. Recent imaging with the Spitzer 
Space Telescope (Hines et al. 2004) demonstrates 
that MIR emission in Cas A correlates well with 
X-ray and optical emission lines rather than with 
synchrotron regions. These authors attribute the 
MIR energy to thermal emission from dust inside 
the SNR, or associated with emission-line gas in¬ 
side the reverse shock region. This correlation 
may arise from grain formation in regions with 
sufficiently cool gas, militating against a smoothly 
varying distribution of MIR emission. 

Nonthermal emission in the IR and radio would 
lead to a ratio of Ss 43 mhz < 3 (Cohen & 

Green 2001). In fact, the nebular ratio is about 
6 for G313.3-1-00.3. This is closer to the values 
found for PNe, for which the median ratio is 12, 
with a semi-interquartile range of 6, than to Hii 
regions (median 24). If the MIR spectrum were 
that of a PN or an SNR dominated entirely by 
emission lines, as in high-excitation objects, one 
would expect spatial variations across the neb¬ 
ula in the intensity ratios among the IRAC bands 
(implicit in Figure 7). We now investigate the 
integrated intensity ratios of 8.0-/rm/5.8-/rm and 
8.0-^m/3.6-^m in G313.3-1-00.3 and compare these 


with known PNe. If it were a PN in which PAHs 
dominated the spectral energy distribution (SED) 
then band-to-band ratios should be rather uniform 
across the nebula. We eliminated the contribu¬ 
tions from the presumed unrelated, bright point 
source within the ring, and from the central star. 
The nebular mean intrinsic values for these two ra¬ 
tios in G313.3-1-00.3, derived by dereddening the 
observed maps, are 3.1 ± 0.2 and 12 ± 3, and are 
quite uniform over the nebula. These uncertainties 
incorporate the differences in corrections between 
a conventional extinction law and that presented 
by Indebetouw et al. (2005). 

A rough estimate of what to expect for a PN 
with only PAH emission comes from assigning 
the observed 3.6, 5.8, and 8.0-^m PN intensi¬ 
ties purely to the 3.3, 6.2, 7.7, and 8.7-/rm PAH 
bands. The predicted ratios would then be 2.7to'3 
(Band 4/Band 3) and I7I3 (Band 4/Band 1) (Co¬ 
hen et al. 1986: Table 2, for PNe; Cohen et al. 
1989: Tables 3 and 4, for PNe). These are very 
similar to the ratios for G313.3-1-00.3. A more pre¬ 
cise approach is to synthesize IRAC photometry 
from MIR spectra of PNe. Volk & Cohen (1990) 
characterized the 7.7 to 22.7-/rm spectra of 170 
PNe according to the dominance of PAHs, broad 
or narrow dust emission features, and of atomic 
fine structure lines of low, high, or very high ex¬ 
citation. The recognition that, in some bands, al¬ 
most all the MIR emission in a PN can be con¬ 
tributed by pure rotational lines of H2 (e.g. Cox et 
al. 1998) implies that one must additionally con¬ 
sider H2 line emission. To quantify this emission in 
an extreme case of a PN entirely dominated by H2 
lines, we have combined the line data for the Helix 
Nebula (Cox et al. 1998) with those for NGC 6302, 
NGC 7027, and Hb 5 (Table 7 of Bernard-Salas & 
Tielens (2005)). In real PNe, the H2 line emission 
that falls within the IRAC bands is also accom¬ 
panied by MIR continuum, fine structure lines, 
and/or PAHs. We have integrated the new (July 
2004) IRAC relative spectral response curves over 
the spectra of 63 PNe observed with the ISO Short 
Wavelength Spectrometer (including three of the 
H2-emitting PNe cited above), and have formed 
the observed intensity ratios 8.0-/im/5.8-/im and 
8.0-/im/3.6-/im. Allowance for extinction in these 
PNe would change the observed ratios very little. 
(For a PN with a modest Ay of d'", the intrinsic 
ratios would differ from those observed by at most 
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4% percent, whichever reddening law is used. In 
fact, Band 4/Band 3 is unaltered with the Indebe- 
touw et al. (2005) extinction.) We have averaged 
the ratios observed for PNe grouped according to 
their most conspicuous MIR spectral character. 

For 16 PNe dominated by PAH emission the 
mean ratios of 8.0-/im/5.8-/xm and 8.0-/xm/3.6- 
/im are found to be 2.9±0.2 and 14±3, entirely 
consistent (within Ict) with those of G313.3+00.3. 
Could another emission mechanism mimic these 
radiance ratios? The most restrictive criterion is 
the 8.0-/im/3.6-/xm ratio. Only for PNe dominated 
by very high excitation emission lines does this 
ratio attain a value (14) comparable to that of 
G313.3+00.3 but, for such PNe, 8.0-/rm/5.8-/rm is 
much too low ('^ 1). No single emission process 
other than PAHs can simultaneously match the 
intrinsic Band 4/Band 3 and Band 4/Band 1 ra¬ 
tios of G313.3-1-00.3. Specihcally, H2 line emission 
alone gives mean ratios of 5.5±1.2 and 4.9±1.8, 
respectively, quite inconsistent with G313.3-1-00.3. 
Might a combination of emission processes match 
the ratios in G313.3-1-00.3? Adding the pure H2 
to the very high excitation line spectrum in the 
proportion of 1:3 would formally produce spatially 
integrated nebular ratios of 2.8±0.5 and 11.9±2.8. 
However, it is impossible for such values to be rel¬ 
atively uniform across a PN because each process 
requires radically different physical conditions. H2 
emission arises near the surface of PDRs but the 
chief cooling lines are in the far-IR, from low exci¬ 
tation ions such as [Cii], [Oi], and [Sin] (Bernard- 
Salas & Tielens 2005), not the species that char¬ 
acterize IRAC’s wavelength range in very high ex¬ 
citation PNe such as NGG 2440 (Bernard-Salas et 
al. 2002), e.g. [Arvi] and [Nevi]. The band-to- 
band ratios of G313.3-1-00.3 are, therefore, consis¬ 
tent only with a PAH-dominated spectrum. 

Figure 8 shows the spatially-integrated SED for 
the nebula, both observed, and after dereddening 
corresponding to an Ay of 36™, using Indebetouw 
et al. (2005). For each data point, asymmetric 
horizontal bars represent the extent of the FWHM 
of the corresponding filter band. The obvious fea¬ 
tures appear to correspond with PAH emission 
bands at 3.3 and 7.7 /rm. Only a weak PAH con¬ 
tribution would be expected in the 14.6-/im MSX 
band from G-G-G modes near 16 /rm. The domi¬ 
nance of the 3.3-/im PAH band in this plot of spe¬ 
cific intensity is misleading. Both this feature and 


the IRAC 3.6-/im bandwidth are far narrower than 
the 7.7 fim band and the 8.0-/im filter, greatly re¬ 
ducing the integrated in-band flux contributed by 
PAHs to the IRAG Band-1 data point as compared 
with Band-4. Between 3 and 30 /rm the ring emits 
only about 700 Lq, even after dereddening by 36™. 
However, the complete absence of any optical or 
NIR detection of the ring renders any meaningful 
estimate of its bolometric luminosity impossible 
from our observations. The bulk of the expected 
luminosity should occur in the optical and UV, 
but we lack any data to deredden at these short 
wavelengths. 

Scaling down the predicted IR thermal emission 
of the entire nebula to that of the central source 
by the ratio of their 13-cm flux densities suggests 
approximately 0.2 and 0.4 mJy in IRAG Bands 1 
and 2. These predictions are of the same order as 
the values measured for the central star, and the 
predicted flux density ratio is in accord with that 
observed, at the la level. (In the above simplistic 
calculation we have neglected the probable higher 
values of Te and ionic charge expected within the 
stellar wind as opposed to the overall nebula.) It 
is somewhat unusual for a star to have a radio 
counterpart. But this a plausible explanation in 
the context of the central star of a PN. Taylor et 
al. (1987) modeled thermal radio emission from 
such a stellar wind with a 1 /r^ density profile, and 
matched the observed radio spectra for a sample of 
PNe with slopes of >1. Truncation of these stellar 
winds also steepens radio spectral indices (Marsh 
1975). Therefore, the central star, or at least its 
wind, is detected in both the IR and the radio. 

We conclude that G313.3-1-00.3 is most likely to 
be a PN within heavy obscuration associated with 
our almost tangential view along the Scutum-Grux 
spiral arm. The regular nebular morphology, the 
presence of a radio and IR source at the geometric 
center, the spatially rather uniform ratio-images 
in key IRAC bands, and the quantitative match to 
what is expected for these band ratios from PAH 
emission, all support this classihcation. We antic¬ 
ipate that this PN will be only the Hrst of many 
highly extinguished planetary nebulae to be dis¬ 
covered by GLIMPSE. 
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Fig. 1.— High-resolution 13-cm greyscale image 
of G313.3-1-00.3 overlaid by 20-cm contours at 0.6, 
0.8, 1.0, 1.2, 1.5 mJy beam“^. The position of 
the central star coincides with the white dot sur¬ 
rounded by a small white circle. 
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Fig. 2.— Radio spectral index map of the neb¬ 
ula between 13 and 20 cm, in the range a from 
about —1 to +2. The star’s position is within the 
small white circle. 20-cm contours are as in Fig. 1. 
White pixels represent indeterminate spectral in¬ 
dices. 
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Fig. 3.— Montage of IRAC high-resolution im¬ 
ages of the ring in the 4 bands. Pixel size is 
0.6". Diffraction patterns can be seen around the 
brighter point sources. 
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Fig. 4.— Brightness profiles through the central 
star’s position in the four IRAC high-resolution 
images. Negative distances from the center are to 
the south of the star. 
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Fig. 5.— 13-cm high-resolution greyscale image 
overlaid by white contours of 4.5-^m emission at 
values of 3, 4, 5 MJy sr“^. The lowest contour 
is 4(7 above the background. The star is seen at 
both wavelengths in the center of the clumpy shell 
and is shown as in Fig. 1. Note the good match 
in spatial extent of the outer boundaries in the IR 
and radio maps. 
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Fig. 6.— 13-cm high-resolution greyscale image 
overlaid by white contours of 8.0-/im emission at 
values of 34, 42.5 MJy sr“^. The lowest contour 
is 4(7 above the background, and corresponds to 
34 MJy sr“^, as does the third contour. The sec¬ 
ond contour represents the bright ridgeline of the 
IR ring (42.5 MJy sr“^). Star shown as in Fig. 1. 
Note the greater spatial extent of the IR outer 
boundary compared with the radio emission. 
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Fig. 7.— MSX 8.3-^m greyscale image of Cas A 
overlaid by white contours of 21.3-^m emission at 
values of 0.75, 1.5, 2, 2.5x10-^ W m-^ sr-h The 
lowest contour is 3.5 (t above the background. 
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Fig. 8.— The observed (open squares, dashed line) 
and intrinsic (filled squares, solid line) SEDs for 
the PN before, and after, correction for 36'" of 
visual extinction. Horizontal bars through each 
data point indicate the FWHM of the associated 
filter bands. Note the clear indications of PAH 
emission at 3.6 and 8.0 /rm. 
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Table 1 

Radio and IR flux densities of the ring and central star. 


Waveband 


Flux Density 

The nebula 



MOST 843 MHz [mJy] 


23 ±3 

ATCA 1384 MHz (low-resolution) | 

mJy] 

35 ±5 

ATCA 2496 MHz (low-resolution) | 

mJy] 

25 ±3 

MSX/CB-02 F{8.3^lm) [mJy] 


130 

MSX/CB-03 F(8.3Aim) [mJy] 


125 

MSX/CB-03 F(14.6Aim) [mJy] 


450 

IRAS F{12nm) (Jy) 


< 2.2 

/RA5 F(25Aim) (Jy) 


0.81 ±0.15 

/RA5F(60^m) (Jy) 


< 8.0 

/RA5F(100/xm)(Jy) 


< 140 

IRAC F(3.6/ini [mJy] 


30 

IRAC F(4.5/ini [mJy] 


25 

IRAC F(5.8y^ni [mJy] 


38 

IRAC F(8.0^ni [mJy] 


150 

The Central Star 



ATCA 1384 MHz (high-resolution) 

[mJy] 

2.6 : 

ATCA 2496 MHz (high-resolution) 

[mJy] 

1.0 ±0.05 

IRAC F(3.6/nn [mJy] 


0.30 ±0.06 

IRAC F(4.5/nn [mJy] 


0.43 ±0.08 
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